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Abstract. Most black hole candidate X-ray binaries show Fourier time lags between softer and harder X-rays. The hard 
photons seem to arrive up to a few ms after the soft for a given Fourier frequency of the perturbation. The energy dependence 
of the time lags has a roughly logarithmic behavior. Up to now most theories fail to explain the observed magnitude and Fourier 
frequency dependence of the lags or fail other statistical tests. We show that the time lags can arise from a simple pivoting power 
law model, which creates the logarithmic dependence on the photon energy at once. A pivoting power law arises naturally from 
jet/synchrotron models for the X-ray emission, but may also be applicable to corona models. A hint to the coherence features 
of the light-curves can be obtained from the power spectral density, which can be decomposed into a few broad Lorentzians 
that could arise from a couple of strongly damped oscillators with low quality factors below one. Using small variations of the 
power law index for each Lorentzian separately the lags can be derived analytically. They show the correct Fourier frequency 
dependence of the time lags. If one assumes variations of the power law index by ±0.2 the model can account for the observed 
magnitude of the time lags in Cyg X-l. The model can also be applied to TeV blazars, where a pivoting power law and hard 
lags have been observed directly in some cases. As a further test we calculated the cross- and auto-correlation functions for our 
model, which also show qualitatively the observed behavior. The auto-correlation function for higher energies has a narrower 
peak than at lower energies and the cross-correlation function is asymmetric but peaks nearly at zero. The coherence function 
for the model is in agreement with the observed data in the Fourier regime, where the model is valid. 
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1. Introduction 

The central part of active black holes seems to consist of the 
black hole with an accretion disk surr ounded by a hot corona 
(see e.g., ISunvaev & Trumped ll 9791 or [Haardt & Maraschi 
1991 ) a nd a jet (e.g ., ISpenceJ Il979t iMirabel & Rodriguez! 
1999, or Fender 2001). However, up to now the accretion flow 
of black holes, jets, their connection, and their relative promi- 
nence are not well understood. The most common active black 
holes are active galactic nuclei (AGN) and black hole X-ray 
binaries (BHXRBs). 

To constrain models and physical parameters of these 
objects it is important to access all observable quantities. 
Besides the spectra the variability is of high importance 
as it can reveal information about the central engine and 
its dynamics. Strong variab ility is a common phenomenon! 
for XRBs (see e.g.. Ivan der Knlll989l) . In BHXRBs the X- 
ray emission is commonly explained b y an accretion disk 
and a Com ptonizin g corona (see, e.g. IShap iro et alJ fl976l 
ISunvaev & Trumpet! 1 9791 lHaardt & Mara schill99ll) . but there 
may also be signif icant contributions f ro m synchrotron emis- 
sion from a jet ( Mar koff et alJ l200ll or iFalcke & Biermanr] 
Il999l) . The jet/synchrotron model predicts a rigid power 



law that can only vary in amplitude and in spectral index. 
Variability in Comptonization models can lead to a pow er law 
X-ray spectrum as well (see e.g., Kylafis & Klimis 1987). Here 
we will investigate whether the short term variability of ac- 
tive black holes can be explained with a rigid pivoting power 
law model. We will concentrate on BHXRBs as detailed light- 
curves are available, but applications to AGN are as well pos- 
sible. 

Usually BHXRBs appear in two distinct states: the hard 
state (low flux levels accompanied with a hard power law spec- 
trum) and the soft-state (normally higher flux and a soft X-ray 
spectrum, see e.g.. Ivan der Klisl [l994). In the hard-state a rel- 
ativistic jet ca n usually been seen in radio observations (e.g., 
Fended 1200 4 We will focus our studies on the hard state, 
where the X-ray spectrum is dominated by a power law. 

A BHXRB in the hard state shows significant short time 
(0.1-100 Hz) variability with a roo t mean square (rms) around 
20% (see e.g., Ivan der Klisl fl99l . It is therefore possible 
to make a detailed statistical analysis of the observed light 
curves. The light curves at different photon energies are well 
correlated as the cross-correlation functi on peaks nearly at 
unity. Furthermore, the coherence function ( Vau ghan & Nowakl 
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1 19971) is nearly unity for a wide range of Fourier frequen- 
cies. However, one often observes hard lags, e.g. the hard 
photons lag behin d the soft photons up to a few millisec- 
onds (e.g.. seelNolan et alJl98lllMivamoto & KitamotoT l989. 
iMivamoto et al.Hl 99 ll or lPottschmidt et al.ll2000t for a defini- 
tion of phase lags see below, Eq. 1131 . The existence of hard 
lags has been explained using Comptonization models. Soft 
photons will be repeatedly up-scattered in a large corona, as 
the harder photons need more inverse Compton processes to 
reach their ener gy this results in hard l a gs. For studies using 
coronae see e.g.. | M iyamotq et alj dl99lh.Fj owaketal 1dl999i). 
iMalzac & Jourdainld2000l) . |Poutanenl(l2002l) . or lBottcher etal] 
(2003). As already noted by these authors, this explanation has 
the problem that one needs huge coronae and the Fourier fre- 
quency dependence of the X-ray time lags cannot be repro- 
duced. 

Additionally the observed auto-correlati on is not repro- 
duced well (see e.g., [Maccarone et al .l 120001) . A different ap- 
proach has been made by Ko tov et al J (12001 ). where the authors 
explain the phase lags with the response of the accretion disk 
to perturbations and present a short discussion of the effects of 
a pivoting power law. 

By the term pivoting power law we mean that the X-ray 
spectrum at different times can always be described by a power 
law, which only varies in the power law index and the overall 
intensity. We mostly consider the case where the amplitude and 
the power law index are correlated. 

The idea of a pivoting power law model arises from re- 
cent theoretical and observational results. The spectrum of 
BHXRBs can b e well descri bed using a coupled jet/accretion 
disk model (see Mark off et all200ll) . Here the disk (possibly a 
optically thin accretion d isk, e.g., such as ADAFs and related 
solutions. lNaravan & Yil l995. plus a standard disk) is only vis- 
ible as an additional component in the UV, while the flat spec- 
trum at radio and optical wavelength and the power law in the 
X-rays is created by synchrotron and inverse Compton emis- 
sion from the jet. In particular, the hard X-ray power law is ex- 
plained as optically thin synchrotron emission from a single re- 
gion at a few hundred Schwarzschild radii from the black hole. 
The power law index depends on plasma parameters (e.g., elec- 
tron temperature, adiabatic index), and may therefore respond 
to changes of the jet power and the accretion rate. As the to- 
tal intensity depends on these parameters as well, the flux and 
the power law index should be correlated. The jet/synchrotron 
model therefore suggests that the X-ray emission behaves like 
a pivoting power law. 

Within the jet/disk picture of Mark off et ail d200lh . TeV 
Blaz ars like Mrk 421 or F anaroff-Riley class I radio galax- 
ies |Fanaroff & Ri levlll974 the unbe amed parent p opulation 
of BL Lacs within the unified scheme. Iljrrv & Padovanill 995) 
show many features of BHXRBs in the low/hard state, namely a 
domination of the spectral energy distribution by jet emission. 
The connection of XRBs in the hard state and jet dominated 
AGN is discussed in Falc ke et alJ j2003). Mrk 42 1 , for exam- 
ple, shows ha rd lags and a positive hardness/flux correlation 
(Zha nel200i . The hardness seems to show a hysteresis effect, 
e.g. the power law index seems to respond slightly after the 
variation of the total intensity. If BHXRBs also have a power 



law from their jets, a similar pivoting power law could play an 
important role. Hard lags and positive or negative hardness-flux 
correlati ons have also been fo und in Seyferts and other AGN 
see e.g., Chian g etalJ fcoOO) or lLamer et al l d2003l) . 

A pivoting power law may also be applicable for 
Comptonization models. Analy zing long term varia bility 
(timescales of days) of BHXRBs Zdziars ki et alJ d2003l) sug- 
gest the existence of a pivoting power law with a pivot point 
around 50 ke V and explains the behavior using Comptonization 
in a corona. They find a negative correlation between flux and 
hardness. These long term variations arise probably from a dif- 
ferent source of variabi lity (e.g., the accret ion rate or an other 
unknown parameter, see Hom an et al. 2002) than the short term 
variations studied here (maybe created by mag netohydrody- 
namic instabilities, see IPsaltis & Normanll2002l or other un- 
known sources). Thus, it is yet unclear if such a correlation 
holds for fast variations and the true hard state. 

In this paper we will analyze in a general way the effects 
of a pivoting power law model, where the power law index is 
correlated with the flux. We calculate the effect on the phase 
lags and the auto- and cross-correlation functions, and present 
a Monte Carlo simulati on of the cohe rence function. In addi- 
tion to the work bv lKotov et alJ d200ll) . who also discussed the 
possibility that the power law index is directly correlated with 
the flux, we include a response time for the change of the power 
law index as a function of intensity. 

In Sect. 2 we describe our parameterization and model. 
With these definitions we derive a general analytic solution for 
phase lags and cross-correlation functions for a pivoting power 
law in Sect. 3. In Sect. 4 the analytic result is compared with 
a Monte Carlo simulation. In the last two sections we discuss 
our model in the context of data from Cygnus X-l and present 
our conclusions. 



2. Parameterization of the pivoting power law 
model 

As we try to calculate the time lags with an analytical approx- 
imation it is important to parameterize our pivoting power law 
model around a reference photon energy eo near the observed 
energies. Let the flux S of our source be a function of photon 
energy e and time t 



S(e,t) = A(t) 



(i) 



-a+W 



(1) 



where a represents the constant part of the spectral index while 
/3{t) accounts for the variations. The function A(t) describes the 
flux at the reference energy eo. As we will consider the case 
that A(t) and /3(f) are correlated, the reference energy e will 
not be the pivot point defined by the minimum of the rms. 

If the changes in spectral index are small and we are observ- 
ing photon energies near the reference energy (In (;|)y8(f) « 1) 
we can expand the equation 



S(e,t)=A(f) 



1+0(0 In - 



(2) 
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and find in Fourier space, denoted by S : 



S ( e. cj) - | j \ U{cS) + /3A(cj) Inlj 



(3) 



As we are interested in phase lags, depending on the co- 
herence features of A(t), it is inappropriate to use red noise for 
the light curve. Information on the coherence of the light curve 
can be guessed from the power spectral density (PSD) defined 
as PSD(w) = A*{u>)A{u>), where the star denotes complex con- 
jugation. We note that the PSD of many BHXRBs can be well 
described by a sum of a few broad Lorentzians 



PSD(tu) = 2^,0,(0*), 



with 



AR]Qm 



2 ' 



(4) 



(5) 



where i 


pne Lorentzian can be centered around ai - (see 


iNowakl 


2000, 


Pottschmidt et al. 


I2003L or 


Bellonietal. 2002). 



The quality factor Q is a measure of the full width half max- 
imum (FWHM) Q = 2 7tfwhm • The normalization factor R de- 
scribes the amplitude of the Lorentzian. It is connected to the 
total rms amplitude as 



rms; = Ri 



1 tan-' {-IQi) 

2 n 



(6) 



(see e.g.. IPottschmidt et alJl2003l) . Lorentzians usually arise 
from damped oscillating systems, for example they are used 
to describe the spectral shape of a laser. The quality factor Q 
describes how strong the oscillator is damped, a high Q denotes 
a nearly undamped system with a strongly peaked PSD, while 
a low value for Q yields a highly damped system with an asym- 
metric, weakly peaked PSD. Usually around four Lorentzians 
with a quality factors Q < 1 are needed to fit the PSD of XRBs 
in the low/hard state. As the origin of these broad Lorentzians 
is still unknown we assume that each Lorentzian is created by 
a strongly damped oscillator excited at random times. For ex- 
ample, these oscillators could be due to excitations at different 
locations on the accretion disk with a Fourier frequency defined 
by the Keplerian rotation that may or may not be transfered into 
the jet. Other possible explanations include magnetohydrody- 
namic instabilities (Psaltis & Norman 2002) or jet precession. 

To simplify the discussion we first look at only one broad 
Lorentzian centered around uq. We assume that the variability 
is created by a damped oscillator. To generate our light-curve 
we use a simple shot noise model (for s hot noise models se e 
e.g.. lTerrelllfl972l lLochner et alJll99ll or lNegoro et al]l200ll) . 
Let us first assume that the light curve of the BHXRB can be 
described as F(t), if this oscillator has only been excited at t = 
with a unit excitation. The overall light curve will be a super- 
position of many excitations at random times and amplitudes. 
If A, describes the amplitude of the excitation at the time f, we 
can write 

A(t) =A DC + J] W - n), (7) 

i 

where Adc describes the constant offset of the flux. The ampli- 
tudes Aj and the excitation times f,- are random variables. We 



choose the normalization of F(f) such that {A 2 ) = 1 . Using this 
process we create a light curve that has the observed PSD and 
the coherence properties given by the oscillator. 

It is unclear whether such simple shot noise mo dels can de - 
scribe the light-curves of XRBs (see e.g. JLochner et al Jl99 lh . 
However, we have to disentangle the contributions of the dif- 
ferent broad Lorentzians to apply our model, so it is very hard 
to test a pivoting power law model using observed light curves. 
Shot noise is therefore one of the best possibilities to create 
artificial light curves available. Furthermore, this approach en- 
ables us to give an analytic solution for the phase lags and the 
cross-correlation function. We discuss below the effect of dif- 
ferent coherence properties. The main result is likely to be in- 
dependent of the shot noise assumption. 

Transfered into Fourier-space we find 



A((o) = A DC S(u) + J] AiFi^e 11 



(8) 



The PSD of a complex damped oscillator is a Lorentzian, 
but as we are interested in real solutions for the light curve we 
have to use a linear combination of the real and the imaginary 
part of the damped complex oscillator. The two fundamental 
real solutions are the instantaneously excited oscillator (cosine) 



and 



F c (co) = i (it + £_) 
F s (a>) = ^ (H + - //_) 



the sine combination, where H ± = 



(9) 
(10) 

. is the Fourier 

tL»o-2iy(tJ±Ct>o) 

transform of the complex oscillator with frequency a> - ±u>q. 
We note that the spectral form of the cosine combination F c 
declines with a>~ 2 like the Lorentzian, while the sine term drops 
with or 4 . 

Given the light curve A(t) we have to choose a physical re- 
sponse of the power law index f3(t). Whatever model one uses 
for the X-ray emission (Comptonization or jet model) the spec- 
tral index depends on physical properties near the accreting ob- 
ject. It is therefore likely that /3(f) will respond to changes of the 
accretion rate. As the region of emission has a characteristic 
size, p may not follow A(t) directly, but may respond a bit later. 
It is also possible that the emission mechanism has a response 
time itself (e.g., for the jet model how fast is the particle accel- 
eration mechanism responding). We will take this into account 
by introducing a response time t. The first order approximation 
of p will therefore depend linearly on A(t-r). We assume that t 
will be a small fraction (a) of the period of the center frequency 
of the Lorentzian, i.e. t = oItxIujq with < a < 1. 

The first order approximation form of /3 and the simplest 
form with the described properties is 



Pit) =yA AC (?-T), 



(11) 



where the subscript AC marks the contributions to the light 
curve from the damped oscillators. The constant part of a vari- 
able will be marked with a subscript DC throughout this paper. 
In the case of a positive hardness/flux correlation, as seen in 
Mrk 421 and other Blazars, the parameter y is positive. For a 
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negative correlation, seen in long timescale variablity of some 
BHXRBs, one has to use y < 0. For illustrative purposes we 
first use a positive hardness/flux correlation (y > 0) in our cal- 
culations and discuss the other case in a separate subsection. 
In Fourier representation one finds for /3 

0(to) = y^ T A AC (to)D(to), (12) 

where we included an additional damping factor D(a>). It 
should take into account the damping of the system creating 
the power law, i.e. if to is too big, the oscillations are so fast 
that they average out, and the power law index -a + /3(t) does 
not change anymore. 

3. Analytic results 

3.1. Definition and energy dependence of phase lags 

The phase lag measures the phase change (corresponding to a 
time delay) between the light curves at two photon energies 
in phase space. The Fourier phase lag <p(o>) is defined as the 
argument of the average cross power spectrum 

0(w) = arg<5(ei,w)*5(e 2 )> (13) 

see e.g., iNowak et alJ Jl999l) . The Fourier phase lag can be 
translated to time lags by dividing through to. The time lags 
have a simple interpretation in the time domain, as they mea- 
sure the time difference between an outburst at two different 
photon energies at a given Fourier frequency. To calculate the 
phase lag from a light-curve we start from 



\§(e u a>)\\§(e2,a>)\ ' 



(14) 



Using the linear approximation (Eq.OO the phase lag can be 
evaluated analytically (see also lKotov et al.l20 01 ): 



sine 



3 [[A* + (/Mnn(|)][A + (/M)ln(g)]] 



A*A 



A* A e\ ' 



(15) 



So, if A*(J3A) has a nonzero imaginary part, the phase lag will 
vary with In e as seen in the observations of BHXRBs. The pho- 
ton energy dependence of the phase lag is therefore indepen- 
dent of the response /3(t) of the power law index, but the overall 
magnitude and the pase lag dependence on Fourier frequency 
to depends on the choice of f3. 

3.2. Fourier frequency dependence of the phase lag 

To derive the phase lag we start with Eq. dl5> . In our linear ap- 
proximation the only unknown component is the Fourier trans- 
form of f3A. Using the convolution theorem and the Fourier 
transform of /3 from Eq. ill 2b we get: 

y j e lu ' T A AC (oj')A(to-to')D(to')dto' 
+yA DC A AC (uj)D(to)e mT . (16) 



PA 




d 



0.01 
0.005 

"^\. io. loo. 

Fourier- frequency to / (1/sec) 

Fig. 1. Phase lag dependence on the Fourier frequency. The 
solid line represents the analytic approximation, the dots are 
the results of a Monte Carlo simulation. The deviation of the 
Monte Carlo simulation at lower frequencies is a numerical ef- 
fect. 



For the case D{to) = 1 it is possible to give an analytic 
calculation of the phase lag. The calculation is described in 
Appendix A and we just give the result for an instantaneously 
excited oscillator: 



= yln — (Adc sin(m>) 
et \ 

(A 2 ) 



+ ^3(e'™/ l(w) ) 



(17) 



where f(to) is given in the appendix. The A denote the random 
excitation power of the oscillators. It is important to note that 
the calculations show that one can treat each excitation of the 
oscillator separately - the contributions of different excitations 
average out. 

The first term in Eq. dl7l i is the result of the pivoting power 
law acting on the constant flux, while the second term repre- 
sents the pivoting power law acting on the pulse itself. So the 
second term depends on ^ and will therefore vanish if the A 
are distributed symmetrically around zero. In that case only the 
first term will contribute in that order of the perturbation series, 
but higher orders can be nonzero. 

The behavior of the phase lag is illustrated in Fig.[T](solid 
line). For small tor the phase lags increase linearly with fre- 
quency. This means that the time lags are constant for small 
frequencies. Once tor w |, the phase lags start to oscillate. If t 
is of the order of — this will happen roughly at to > too. 

The result that the phase lag starts to oscillate is at least 
partly due to our deterministic law of a fixed look-back time 
t. A more realistic approach would be that t is itself a ran- 
dom variable, e.g., Gaussian distributed around a mean value. 
As the sine is linear for tor <sc 1 the phase lags for low Fourier 
frequencies will not change. However, for tor 1 the contri- 
butions to the phase lags for different r will average out. The 
oscillation will be further reduced as the PSD of a real system 
is described by several Lorentzians, i.e. for to » a»o a second 
Lorentzian will dominate the first one. It is therefore likely that 
one will never observe the oscillating part of the phase lag. 
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To take the probably statistical nature of r and additional 
damping of the response of the power law into account, we 

choose a damping term D(a>) + 1 , for example D(oj) = e <5 (" , °) , 
which will cut off the phase lags at a given frequency. With an 
appropriately chosen damping the phase lag will not oscillate 
but stay at zero for higher Fourier frequencies. 

Hence, we conclude that the phase lag in the case t > can 
be approximated by 



sin( 



y sin( lot)D((jS) In — 



(18) 



and is independent of the exact shape of the pulses. We find 
hard lags for a positive hardness flux correlation and a positive 
look-back time. If one changes the parameters, e.g., negative 
hardness flux correlation and positive look-back times, soft lags 
can be obtained. 

We can now verify that our simplification made in Eq. @ 
was appropriate. We assumed that all excitations are identical 
and only vary in amplitude and excitation time. In a real system 
each excitation will have a different shape. But as the contribu- 
tions to the phase lags from different excitations average out, 
each pulse contributes as if there are no other excitations (the 
sum in Eq. IA.4I only runs over the diagonal part). The over- 
all phase lag will be the average of all pulses or the phase lag 
of an average shaped pulse. In the zeroth order approximation 
the coherence does not play a role for the phase lags (as long 
as there is a constant flux component). Thus the shot noise as- 
sumption is not a crucial ingredient for the model. The impor- 
tant assumption is only that the PSD can be decomposed into 
different Lorentzians, which do not interact with each other. 

In the case that the power law index responds instantly to a 
change in the accretion rate, the first term of Eq. d!7i will van- 
ish and the phase lag depends to first order in y linearly on (A 3 ). 
The lags would therefore vanish in the case of a symmetric dis- 
tribution of the A (ma ny systems show asymm etric excitations, 
for one example see ISpruit & Kanbachll2002l) . Furthermore, 
they depend strongly on the linear combination used for the 
pulse shape. The lags created by the sine and cosine term are 
shown in Fig. [5] Nevertheless, if the power law index F is vary- 
ing by AF w 0.2, as will be used here later on, the phase lags 
due to the pivoting power law will contribute significantly to 
the observed phase lags. The sign of the lag changes with fre- 
quency as mentioned bv lKotov et alJ (EoOli). who evaluated an 
instantly changing power law for a real light curve. 

3.3. Cross-correlation Function 

The phase lags depend only on the phase of the Fourier- 
transform and contain no information on the amplitude. To 
gain information also on the amplitudes we consider the 
auto-correlation and cross-correlation functions. The cross- 
correlation function is defined as: 

C(e u e2,f) = J" (S(e u t)-(S(eO)) 

(S(e 2 ,f + r)-<S(e 2 )»df. 



(19) 



We have chosen the signs such that the cross-correlation func- 
tion between a lower and a higher photon energy peaks to the 



M 

-J 

s 

Oh 



0.4 
0.3 
0.2 
0.1 
0. 
-0.1 
-0.2 
-0.3 



0.01 0.1 1. 10. 

Fourier-freqency u>/ (1/sec) 



100. 



Fig. 2. Difference between the sine and cosine part. The magni- 
tude of the sine lags (dashed) have been magnified by a factor 
of ten. 

right if we observe hard lags. Expressed in Fourier-space we 
find: 



C(ei,e 2 



5 A c(ei,w) t 5 A c(e2,w)e- lwf dw. (20) 



If we insert the expanded expression Q for S (e, co) and only 
consider terms up to 0(y 2 ) we find: 



C(e u e2,f) = j A AC A AC du> 

+ j In e 2 +/3AA In ei ) Ac e" iwf dw. 



(21) 



For simplicity we have set the reference energy eo = 1 . 

The first component of the integrand represents the auto- 
correlation function of the light curve at the reference energy: 



Ja ac 



A AC e~ lWT doj 



2 P 2 (cJ)eT 



=: P(f). 



(22) 



The second component takes a longer calculation using the 
results for/M of Appendix A and is described in Appendix B. 
With the function S(t, f ) defined in the Appendix we find for 
the cross-correlation function 



C(ei , e 2 , t) = P(f) + y (In e x P(f - t) + In e 2 P(f + t) 
+ In e 2 S(T, f ) + In eiS(T, -f )) . 



(23) 



As described in the Appendix S peaks at f = r and decays 
faster than P{f). 

Whether the auto-correlation function has a steeper peak 
for higher photon energies e depends on the parameters. While 
the terms in the first line of Eq. ( I23> makes the auto-correlation 
function broader for higher energies (~ log e) the terms in the 
second line have the opposite effect. For small Q and t the peak 
is steeper for higher photon energies while for larger values of 
t and Q the opposite effect is found. 

If the excitations of the oscillator are symmetric, i.e. 
(A 3 ) = 0, S vanishes, and only higher order terms contribute 
to the cross-correlation function. The feature that the cross- 
correlation function can have a steeper component for higher 
energies remains in this case as seen in the Monte Carlo simu- 
lation. 
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Fig. 3. Numerical evaluation of the calculated cross-correlation 
function for one Lorentzian only, solid line: low energy auto- 
correlation function, dashed: higher energy ; dotted: Cross- 
correlation function 
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Fig. 4. Monte Carlo simulation of the calculated cross- 
correlation function for a negative hardness/flux correlation and 
a negative look-back time. 



A numerical evaluation of the cross correlation function 
with a damping factor of unity is shown in Fig. Q- The cal- 
culated auto-correlation function is, for the parameters used, 
steeper for higher energy photons. The cross-correlation func- 
tion between two energies lies in between and is slightly asym- 
metric. It is important to note that whether the auto-correlation 
function is steeper for higher photon energies depends strongly 
on the parameters, e.g. the look-back time t. 

If one includes a damping factor D{cS) that damps higher 
frequencies in the response of the power law index, the light 
curve itself will have less power in high Fourier frequencies for 
higher energies. Therefore the auto-correlation function will 
have a flatter peak for higher energies than without a damp- 
ing factor. The effect described above (steeper auto-correlation 
function for higher energies due to the >S-term) and the effect 
of the damping factor can cancel each other partly. 

3.4. Negative correlation of hardness and flux 

In the previous sections we have considered a pivoting power 
law model with a positive hardness flux correlat ion (y > 0) and 
a positive look-back time t. However IZdziarski et alJ (|2003) 
observed a negative hardness flux correlation and a pivot point 
around 50 keV for Cyg X-l in long term variablity (timescale 
of day s). On shorter timescales lLi et alJ i 19991) and lFeng et alJ 
(1999) report a negative correlation as well. Therefore a nega- 
tive correlation can also be present on the short timescales dis- 
cussed here. A negative hardness/flux correlation (y < 0) and a 
positive look-back time (r > 0) leads to soft lags. Soft lags have 
for ex ample been observed in X-ray burst oscillations (Ford 
119991) . However, the observed lags are approximately 1 rad, 
much larger than the lags discussed here. So this effect may 
be due to other mechanisms. 

In order to create hard lags with a negative hardness/flux 
correlation (y < 0) one has to use negative look-back times r. 
This means that the power law index is changed slightly before 
the flux changes. Now the pivot point defined by the rms is at 
higher photon energies as the reference frequency eg. If eo « 3 
keV pivot points around 50 keV can be reached. In this case 



two signs are changed in Eq. Jl 8i resulting in hard lags as in 
the case of y > 0. If one changes the sign of r we have to 
consider that the analytic behavior of Eq. d!6i changes as well. 
The first term of the phase lag (sin ~ y sin(wr)) is the same 
for negative and positive correlation. However, the second term 
does not depend on 3e"" T as before (Eq. IA.9l see Appendix A). 
We find instead 



sine 



y in | Ja dc sin(rw) + ^ 3 (f 2 (co))\ , (24) 



where fo is a rational function in o>. If the constant flux Adc 
dominates the overall flux the non oscillating second term will 
only be a minor correction. The phase lag from a pivoting 
power law with negative look-back times t and negative cor- 
relation will therefore yield similar phase lags as before. 

Even though the second term of the phase lags has a dif- 
ferent structure, the auto-correlation function has the same an- 
alytic structure as in the case t > (Eq. 1231 . Only the func- 
tion S has to be changed. As y is now negative, the first line 
of Eq. d23t already yields a steeper auto-correlation function 
for higher photon energies e. As the second line depends on 
y as well, it can broaden the auto-correlation function. Which 
one of these effects dominates depends on the parameters used. 
Compared to the case with a positive hardness flux correlation, 
the auto-correlation function has a steeper decline for small val- 
ues of f (see Fig.|4} and a slower one for larger f. 

In summary, the two cases positive hardness/flux correla- 
tion and positive look-back time and the case with a negative 
hardness/flux correlation and a negative look-back time have 
similar hard phase lags. Only the parameters have to be ad- 
justed somewhat differently. However, while in the first case 
the auto-correlation function is flatter for small lags and steep- 
ens for larger lags, the opposite effect happens in the other case. 
With appropriately chosen parameters both possibilities seem 
to be able to reproduce the qualitative behavior seen in the ob- 
servations of Cyg X-l by Ma ccarone et al.l (12000). 
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4. Monte Carlo Simulations 

The analytic results have used a linear approximation for the 
response of the power law. For observations near the reference 
photon energy — 1 and small responses of the power law 
index y« 1 this approximation will be valid. However, as one 
often observes photon energies with Sffi > 10 we created a 
Monte Carlo simulation. This will also enable us to consider a 
system with more than one broad Lorentzian. The light curve 
for one Lorentzian was generated as described by Eq. @. The 
random variables A, are chosen to be the absolute value of a 
Gaussian distribution and the excitation times f; are uniformly 
distributed. Each individual pulse is described by the sum of 
the cosine and sine damped oscillator (F = F c + F s ). Once 
the light curve at the reference photon energy is created we can 
calculate light curves at other photon energies using Eq. Q and 
derive the phase lags. The comparison between the numerical 
results and the analytic calculation is shown in Fig.^ We find 
that the analytic first order approximation and the Monte Carlo 
simulation are in good agreement. 

If one assumes that the different oscillators creating the 
broad Lorentzians do not interact with each other one can 
calculate a light curve for a system described by multiple 
Lorentzians by creating light curves for each Lorentzian sep- 
arately and then superpose these light curves to get the over- 
all light curve. This is used in the application to Cyg X-l, in 
Sect.|5] 

Besides the phase lags we can calculate the cross- and auto- 
correlation functions of the light curve. The Monte Carlo sim- 
ulation reproduces the effect seen in the analytic calculations: 
for the parameters used the auto-correlation function is steeper 
for higher energy photons than for the low energetic once. In 
Fig. 10} we show the result of the Monte Carlo simulation with 
a negative look-back time. The cross-correlation function is 
asymmetric as expected for hard lags. 



E 
P 



Cyg X-l time lag (2.5 Hz) 
Pivoting power law 




10 

log Photon Energy / keV 

Fig. 5. Time Lags versus photon energy at a Fourier frequency 
of 2.5 Hz. The solid curve represents the analytic solution for 
the pivoting power law model with a reference energy of 2.0 
keV. 




Fig. 6. Illustration for the superposition of the different broad 
Lorentzians. In the superposition every Lorentzian should dom- 
inate the phase lag around its center frequency. 



5. Applications to BHXRBs 

We are now able to apply the model to the BHXRB Cygnus X- 
1 . Throughout t his section we u se d ata for Cyg X- 1 wh ich was 
taken from lNowak etail J 19991) and lKotov et all fcOOll) for the 
photon energy dependence of the phase lag. 



5.1. Photon energy dependence of the phase lags 

The pivoting power law model predicts that the energy de- 
pendence of the phase lags is logarithmic (Eq. II 51 . In Fig. [5] 
we show the phase lags measured with a constant Fourier fre- 
quency of 2.5 Hz (dots). The solid line represents the pivoting 
power law model. The reference energy e used for the fit is 
2.0 keV. For lower photon energies the accretion disk will be- 
come increasingly important and will start to dominate over the 
power law component. In this regime a simple pivoting power 
law will not be sufficient to describe the behavior and devia- 
tions from the logarithmic energy dependence of the lags are 
likely. 



5.2. Fourier frequency dependence of the phase lags 

The PSD of Cyg X- l can be wel l fitted with four b road 



Lorentzian s (see e.g.. 



IPottschmidt et alJl2003 



Nowak 2000, Belloni et al. 2002 and 



J. The phase lag from each Lorentzian 
will dominate the overall phase lag around its center fre- 
quency. The idea is illustrated in Fig. |6] where we calculated 
the phase lags for each Lorentzian separately. The look-back 
time r is always the same fraction of the center frequency 
of the Lorentzian, which reduces the number of free param- 
eters. In this plot we assumed that the excitations are similar 
to those of an instantaneously excited oscillator and used as 
damping factor a Lorentzian with Q — 0.25 (Center frequen- 
cies / = 0.25, 1.2, 7.3 Hz, t = 0.03/- 1 ). 

If we assume that the different broad Lorentzians do not in- 
teract with each other we can do the superposition with our 
Monte Carlo simulation. First we generate a separate light 
curve for each Lorentzian. In this way we have the possibility 
to consider that each broad Lorentzian has a different look-back 
time r, which we assumed to be a fixed fraction of the period 
of that Lorentzian. The overall light curve is the superposition 
of the four light curves. The PSD of the artificial light curve is 



1 

0.1 
0.01 
0.001 
0.0001 
le-05 
le-06 
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Table 1. The parameters used for Cyg X-l 




artificial light curve 

Lorentzian 1 

Lorentzian 2 

Lorentzian 3 

Lorentzian 4 

Sum 

■ ■' , i_J , i_ 



0.001 



0.01 



0.1 

Frequency / Hz 



Fig. 7. The PSD of the artificial light curve together with the 
four broad Lorentzians. For their parameters see Tabled 



shown in Fig. The parameters of the Lorentzians are given 
in Table [2 The pulse of the shot noise is a superposition of the 
damped sine and cosine oscillator, so its PSD is not exactly a 
Lorentzian, but this only leads to small deviations of the overall 
PSD. The PSD has been used to find the overall normalization 
of the shot noise. 

The damping factor D(a>) used in the Monte Carlo simu- 
lation is a Lorentzian centered around the center-frequency of 
oscillator creating the broad Lorentzians with the half quality 
factor as the oscillator itself. The constant power law index a 
has been fixed to 1.7, and the correlation factor y = 0.26. The 
look-back time was chosen to be r = 0.04/ -1 , where / denotes 
the center frequency of the Lorentzian. In Fig.[8]the phase lags 
from this Monte Carlo simu lation are plotted i n comparison 
with the observed values from Nowak et al. (1999). The param- 
eters are summarized in Tabled The peak frequencies, where 
the Lorentzians contribute mostly to the rms, are also given in 
this table. The peak frequency (/ pea k — f Vl/(4Q 2 ) + 1) is of- 
ten used for frequency correlations (see e.g., Nowak 2000). The 
parameters of the Lorentzians are wit hin the range of the val- 
ues given by Pottschmi dt et alJ (12003). The effect of the fourth 
Lorentzian on the phase lags is rather small, as this Lorentzian 
does not dominate the overall PSD at its center frequency. 

The used parameters are not unique. For example, a small 
change of the hardness/flux correlation parameter y can be 
compensated by a change of the look-back time t. Furthermore, 
the central frequencies of the Lorentzians needed to fit the 
phase lags depend on the damping factors and the positions of 
the other Lorentzians. These parameters where chosen to give a 
steeper auto-correlation for higher photon energies (see below). 
If one would choose different fractions of the center frequency 
of the Lorentzian for the look-back time r and allows for a dif- 
ferent hardness/ flux correlation factor y for each Lorentzian 
better fits would be possible. 

A similar result for the phase lags can be obtained us- 
ing a negative hardness flux correlation and a negative look- 
back time. One possibility would be to choose y = -0.16 and 
r = — 0.06/ -1 and to center the Lorentzians given in Table^at 
slightly different positions (/ = 0.15, 1.2, 6.0,25Hz). 
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Fig. 8. Phase-lag dependence (0-3.9keV to 8.2-14.1 keV) on 
frequency. Four Lorentzians with parameters given in Tabled 
have been supe rposed with our Mo nte Carlo Code. The data 



nave been superposed witn our M 
were taken from lNowak et alJ (1 199' 



5.3. Auto- and cross-correlation function 

In addition to the phase lags we evaluated the auto- and cross- 
correlation function for the simulated light curves at differ- 
ent photon energies. In Fig. [9] we show the auto- and cross- 
correlation function for the parameters used in Fig. [8] here 
with a positive hardness/flux correlation and a positive look- 
back time. With these parameters the auto-correlation function 
is steeper for higher photon energies and the cross-correlation 
is slightly asymmetric. The plot shows qualitat ively the same 
behavi or as the observed correlation shown in Macca rone et alJ 
(2000). However, this behavior depends strongly on the param- 
eters, such as the look-back time t or even the relative promi- 
nence of the different Lorentzians. 



5.4. Coherence function 

The coherence function measures the linear correlation of the 
two light curves at different photon energi es. For an introduc- 
tion a nd a discussion of its properties see Vau ghan & Nowakl 
(1997). It is defined as 



2 \(S*(e u aj)S(e2,aj)) d \< 



(25) 



where denotes an average over different realizations of 
the statistical process, e.g. different light-curves. If the coher- 
ence function is unity the light-curves at different photon en- 
ergies are connected by a linear transformation. To evaluate 
the coherence function for the pivoting power law model we 
simulated 100 different 5000 second light curves with a time 



Elmar Kording and Heino Falcke: X-ray Variability of BHXRBs 



9 



0.1 



0.01 



0.001 



11 keV 
2keV 
Cross-correlation 



_i_ 



_i_ 



_i_ 



_i_ 



-1.5 -1 -0.5 0.5 1 
Time lag / sec 



1.5 



Fig. 9. Auto- and cross-correlation function created by the 
Monte Carlo simulation. See the Tabled for the used param- 
eters. 
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Fig. 11. Comparison of the phase lags from a normal set of 
Lorentzians and a 'failed state transition', where the Lorentzian 
two dominates. 
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Fig. 10. Coherence Function for a pivoting power law model. 
See the Table[T]for the used parameters. 



resolution of 5 ms. The result is shown in Fig. [10] together 
with the measured coherence function of Cyg X- 1 . The high- 
est frequency Lorentzian included in our calculation is centered 
around 25 Hz. Due to our damping factors the two light-curves 
are therefore perfectly correlated at higher Fourier frequencies. 
This leads to the rise of (10) for higher Fourier frequencies to 
one. The rise is therefore an artifact of our choice of the damp- 
ing factors and probably missing higher Lorentzians. In the low 
Fourier frequency domain (< 0.01 Hz) the PSD i s usual ly dom- 
inated by a power law noise component (Nowak 2000), which 
may be responsible for the lower coherence function in this 
regime. We conclude that the model is consistent with the data 
in the Fourier frequency regime where the model is valid. 

5.5. Failed state transitions 

Pottschmidt et al. (2000 & 2002) identified some flares of 
Cyg X-l as 'failed state transitions'. During these flares 
the X-ray spectrum softens and that the rms amplitude de- 
creases. The PSD is then usually dominated only by the second 
Lorentzian (peak frequency m 3 Hz) and to some extend the 
third Lorentzian (peak frequency m 9 Hz). The phase lags dur- 



ing the 'failed state transitions' increases significantly in the 3- 
10 Hz regime, while it stays nearly constant compared to a nor- 
mal hard state below 1 Hz and above 10 Hz. The third and sec- 
ond Lorentzian dominate the PSD between 3 and 10 Hz. This 
leads these authors to the conclusion that the increased phase 
lag can be attributed to these Lorentzi ans. It has already been 
pointed out bv lNowakl feOOOh and lPottschmidt et al.l d2003) that 
the phase lags could be reduced due to the superposition of the 
different Lorentzians. If one Lorentzian starts to dominate the 
overall PSD the phase lags will therefore increase. For exam- 
ple in Fig. |n]we reduced the strength of the first and third 
Lorentzian (R = 0.1) while leaving all other parameters un- 
changed. This leads to an increase of the phase lags, where 
the second Lorentzian dominates, but decreases the lags fur- 
ther away. The increase and decrease of the phase lags depend 
on the other model parameters as well. While inc reasing lags in 
the 3-10 Hz regime are observed JPottschmidt et al.l2003l) . the 
decrease can not be found in their plot. However, during state 
transitions other model parameters may change. It may be that 
the look-back time t or the hardness flux correlation increases, 
resulting in larger lags. To fully understand this intriguing phe- 
nomenon a more detailed study is needed. 

6. Summary and Conclusions 

In this paper we have discussed the effect of a variable pivot- 
ing power law in the spectrum of an astrophysical source on 
its timing beh avior, particularly fo r the Fourier phase lags (see 
e.g., [Miyamot o & Kitamotolll989i) . the cross/auto-correlation 
and coherence function. This model is applied to black hole 
X-ray binaries. From this approach follows immediately that 
the Fourier p hase lag dependen ce on photon energy is logarith- 
mic (see also Kotov et al. l200lh . which is observed in Cyg X- 1 . 
This result is independent of the choice for the response of the 
power law and the coherence features of the light curve. 

To derive the Fourier frequency dependence of the phase 
lags, the coherence properties of the light curve are needed. 

Hints to the coherence of XRBs can be foun d in the PSD, 

I it -i 

which can be fitted by a few broad Lorentzians (Nowak 2000, 
Pottschmidt et al. 2003 or Bello ni et al.l 120021) . A Lorentzian 
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normally arises from a damped oscillator. Therefore, we as- 
sume that the variations of the light-curve are generated by 
randomly excited damped oscilla tors, i.e. we use a simple shot 
noise model (see e.g. jLochner et alJl99lh to generate our light 
curve. The power law index was chosen to respond linearly to 
flux changes including a response time (t), i.e. the power law 
index responds slightly after or before total intensity changes. 

The analytic calculations reduce the Fourier phase lags to 
a simple expression, sin <p ~ y sin(wT) In ~ for one Lorentzian, 
where y is the flux/hardness correlation parameter and a> is the 
Fourier frequency. This law will break down around the center 
frequency of the Lorentzian, due to an included damping of the 
response of the power law and stochastic variations of the look- 
back time t. The phase lag will therfore simply drop to zero for 
Fourier frequencies much higher than the center frequency of 
the Lorentzian. 

To obtain hard lags one has to use a positive hardness/flux 
correlation and a positive response time (response after the 
change of the flux) or a negative hardness/flux correlation and 
a negative response time. If the power law index changes by 
Ar w 0.2 around F = 1 .7 and the look-back time is of the order 
of 10% of the period of the excitation we can account for the 
observed magnitude of the phase lags in Cyg X- 1 . Soft lags can 
be achieved by using a negative hardness/flux correlation and a 
positive response time - or vice versa. 

The result for the phase lags of one Lorentzian is fairly in- 
dependent of the exact shape of the excitations. However, the 
coherence properties of the light curve become more important 
if one superposes different Lorentzians, as needed for Cyg X- 1 . 
Using four Lorentzians we were able to rep roduce the observed 
hard lags of Cyg X-l tNowak et alJI 19991) using parameters of 
the Lorentzians within the published range. Similar hard lags 
are observed for BL Lacs (Zhang 120021) . where the pivoting 
power law model may be applicable as well. The superposi- 
tion of the Lorentzians is likely to p lay an important role i n 
the 'failed state transitions' found by Pottschmi dt et alJ 112000). 
During these events one Lorentzians normally dominates the 
overall PSD, the effect of the superposition is reduced resulting 
in larger lags. 

If one does not allow for a look-back time, a pivoting power 
law will nevertheless create phase lags of the order of magni- 
tude of the observed values. However, these phase lags change 
their sign with Fourier frequency, which is not seen in the data. 
However, when using a pivoting power law model as some- 
times used to explain the rms behavior (see e.g., Zd ziarski et aD 
120031) . one has to take these lags into account. 

Besides the phase lags we also calculate d the auto- an d 
cross-correlation function (see e.g.. iMaccarone et alJ EoOO). 
They show the qualitative correct behavior seen for Cyg X- 1 . 
It is important to note, that while the result for the phase lags is 
fairly independent of the parameters and the form of the damp- 
ing factor, the auto-correlation function can change its qualita- 
tive behavior. For example, if one uses large look-back times 
and large quality factors Q the higher energy auto-correlation 
function will have a broader peak than the one for lower energy, 
the opposite of what observations suggest forXRBs. 

The coherence function, which measures the linear cor- 
relation between the light-curves at different photon energies 



JVaughan & Nowa3l 19971) . has been calculated and compared 
with observations of Cyg X-l. The model is in agreement with 
the observations, for Fourier frequencies where the numerical 
model is valid. 

Both models (the jet/synchrotron model and the 
disk/corona models) predict the existence of a pivoting 
power law. At least the jet/synchrotron model creates a rigid 
power law without spectral breaks. We have shown that such a 
rigid power law model is consistent with the data of Cyg X-l. 
A more detailed analysis using probably an original light-curve 
and physical parameters for both models is needed. We also 
point out, that Cyg X- 1 cannot be de scribed by a simple power 
law (see e.g.. iDi Salvo et al]l200ll) . indicating contributions 
from different emission regimes that most certainly will also 
complicate the timing behavior. 

We conclude that with a rather simple ansatz for a pivoting 
power law model we can explain many of the complex fea- 
tures in the phase lags, cross- correlation function and coher- 
ence function seen in the hard power law emission of XRBs. 

Acknowledgements. The authors would like to thank an anonymous 
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Appendix A: Derivation of the phase lags 

The Fourier transform of a damped oscillator is a sum of two 
similar terms. To simplify the calculation let us first consider 
only the first term 

2/?V3^o~ 



(A.l) 



(jJq - 2iQ(a> - a>o) 

In order to evaluate the equation for the phase lags d!5l > we have 
to calculate the convolution needed for/M (Eq.I16>: 

/3A = yj ^F^F^co-a}')^^^'''^'^^' 

+ A DC F 1 (w)J]V " 

(A.2) 

The excitation times f, are random variables, therefore the sum 
over the non-diagonal elements will yield zero as the random 
phases average out. 

The integral can be solved using complex analysis. It has 
two poles: one from F i (u') at w' = ojq - i |^ and the other from 
F\(lo - a)') near u>' — u - cjq + i|^. If t is positive, we could 
close the integral path with a half circle for positive imaginary 
parts. This path only includes the pole at a>' = co-coo+i^. Now 
we can apply the residual theorem, and find for the integral 

2nRloJa _ e -i+-(— o). (A3) 



u>o - iQ(to - 2u)q) 



Inserting this expression in Eq. dl5> we get 



sin (f> — y In — ( Adc sin(m>)+ 



<j / 2nR 2 k> a e -T^+iT(fc^t^o) p 
y 6j -i2(w-2a) ) 



(A.4) 



The real instantaneously excited, damped oscillator is a 
sum of two terms like F\{to). The result will therefore consist 
of four terms, which can be calculated in analogy. We find for 
a positive look-back time t: 



sin d> — yln — (Adc sin(m>)+ 



with 



fm = 



nR 2 aj Q F*(aj)e' T ^ 



(A.5) 



2F*F(a>) \tJo - iQ(w - 2ll>q) 
gi™ (e'™ + e _i ™°)\ 



(A.6) 



— i2(<y + 2u)q) coo — iQco 



The other real solution for the damped oscillator is the sine 
term (see Eq. For this solution one has to change the two 
last signs and include a factor -i. The only fast changing part 
of the second term of Eq. iA.5\ is e 1 ™. Therefore, the dominant 
part will be sin(rw) as the first term. 

In the second case, t < we have to close the integral path 
in Eq. iA.2\ with a half circle in the negative imaginary plane. 
The integral path now includes the pole at a>' = coo - i|g and 
the residual theorem yields: 



2nR 2 u)o 



ojq - iQ(o) - 2o) ) 
and we find for the phase lags: 



(A.7) 



sin0 = yln — Ape sin(Tw)+ 



y wo-ig(oi-2wo) iv ' i J 



(A.8) 



The real instantaneously excited oscillator yields for this 
case: 

sin0 = yln | (a dc sin(™) + ^3 (/ 2 (w))J , (A.9) 



with 



2F*F(u>) \u>o- iQ(w - 2a>o) 
e -"^F(o)) (e- 1 ™ + e i ™»)F(w) ^ 
wo - iQ(d> + 2a>o) too - iQw I 



(A. 10) 
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Appendix B: Derivation of the cross-correlation 
function 

In order to evaluate the second integral of the cross-correlation 
function (Eq. 1211 

j (A*p~A In ei+p~AA In ei) AC e - ^ dw, (B.l) 

we start with Eq. \A.2\ for f)A. Using this result it is possible to 
evaluate the integral with the same arguments as used for f}A. 
We define three supplementary functions 



g-irwo g-i(2T-T)«o 

<Hi(t,t) = — — + 



e 



i(T-T)O) 



• 2Cos(tu)q)- — — + cc. 

iojQ - 2iQo>q 

giro>o-2ifaJo gir6Jo-2ifwo 



3u>o + 2iQu>o 3u>o + 6iQa>o 
1 

• 2Cos(t(±>o)- — — + cc, 

3w - 2igw 



(B.2) 



(B.3) 



where cc denotes the complex conjugate, and 

S(t,t) = fl^'We^ (© T _ f e- ( - f) i 
Wi(t, f) + © f _ r e 2(T - f) i<H 2 (T, f )) J] <l 



(B.4) 



where v is the step function. It is defined as A = for x < 
and © A = 1 otherwise. We find for the cross-correlation func- 
tion: 

C(ei , e 2 , t) = P(f) + y(ln £iP(f - t) + In (f + r))+ 

y (In e 2 S (t, f) + In ei S (r, -f )) . ' B ' 

The first line of this equation will broaden the auto-correlation 
function (ej = e£) for higher photon energies, positive look- 
back times r, and positive hardness flux correlation y > 0. 
However, S{t, f) peaks at f = r and will decline faster for 
higher r > f than P(r) (falls with e~ T ) due to the factor e 2(T ~ T \ 
This leads to steeper auto-correlation functions for higher e. 
Which effect dominates depends on the parameters used. 

For negative look-back times t and y < we have to use 
Eq. JA.7t for the integral of Eq. JA.2> resulting in a different 
cross-correlation function. The overall form will stay the same, 
as Eq. JA.7t differs from Eq. JA.3> only by the argument in 
the exponential function - the analytic behavior (poles) is the 
same. Here the first line of Eq. JB.5I leads to a steeper auto- 
correlation function (y < 0) while the second line has the op- 
posite effect. 



